Astronomy & Astrophysics manuscript no. TTauS' variability accepted 


©ESO 2010 


April 21, 2010 





Variable accretion as a mechanism for brightness variations in 

TTauS* 

R. van Boekel\ A. Juhasz\ Th. Henning^ R. Kohler^, Th. Ratzka^ "*, T. Herbst^'^, J. Bouwman^, and W. Kley^ 
' Max-Planck Institut fiir Astronomie, Konigstuhl 17, D-69117 Heidelberg, Germany 

- Landessternwarte, Zentrum fiir Astronomie der Universitat Heidelberg, Konigstuhl 12, D-691 17 Heidelberg, Germany 
^ Universitats-Sternwarte Miinchen, ScheinerstraBe 1, D-81679 Miinchen, Germany 

* Astrophysikalisches Institut Potsdam, An der Stern warte 16, D- 14482 Potsdam, Germany 

^ Herzberg Institute of Astrophysics, National Research Council of Canada, 5071 West Saanich Rd, Victoria, BC, V9E 2E7 Canada 

* Institut fiir Astronomie & Astrophysik, Universitat Tiibingen, Auf der Morgenstelle 10, D-72076 Tiibingen, Germany 

Received 2009 December ; accepted 

ABSTRACT 

Context. The southern "infrared companion" of TTau is known to show strong photometric variations of several magnitudes on 
timescales of years, as well as more modest < 1 mag variations on timescales as short as one week. The physical mechanism driving 
these variations is debated, intrinsic luminosity variations due to a variable accretion rate were initially proposed, but later challenged 
in favor of apparent fluctuations due to time- variable foreground extinction. 

Aims. We seek to investigate the nature of the observed photometric variability. Based on simple geometric arguments and basic 
physics laws, a minimum variability timescale can be derived for which variable extinction is a viable mechanism. Because this 
timescale increases rapidly with wavelength, observations at long wavelengths provide the strongest constraints. 
Methods. We used VISIR at the VLT to image the TTau system at two epochs in February 2008, separated by 3.94 days. In addition 
we compiled an extensive set of near- and mid-infrared photometric data from the literature, supplemented by a number of previously 
unpublished measurements, and constructed light curves for the various system components. We constructed a 2D radiative transfer 
model for the disk of T Tau Sa, consisting of a passively irradiated dusty outer part and a central, actively accreting component. 
Results. Our VISIR data reveal a +26±2% change in the TTauS/TTauN flux ratio at 12.8/im within four days, which can be 
attributed to a brightening of T Tau Sa. Variable extinction can be excluded as a viable mechanism for the observed flux variation 
based on the short timescale and the long observing wavelength. We show that also the large long-term photometric variability and its 
associated color-magnitude behavior can be plausibly explained with variable accretion. However, variable extinction is also a viable 
mechanism for the long-term variability, and a combination of both mechanisms may be required to explain the collective photometric 
variability of Sa. 

Conclusions. We conclude that the observed short-term variability is caused by a variable accretion luminosity in T Tau Sa, which 
leads to substantial fluctuations in the irradiation of the disk surface and thus induces rapid variations in the disk surface temperature 
and IR brightness. Both variable accretion and variable foreground extinction can plausibly explain the long-term color and brightness 
variations. We suggest that the periods of high and variable brightness of Sa that we witnessed in the early and late 1990s were due to 
enhanced accretion induced by the periastron passage of Sb, which gravitationally perturbed the Sa disk. 
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1. Introduction 



Six decades ago Ambartsumian| (1947 1949[ l first proposed 
T Tauri stars to be newly formed low-mass stars. They had been 
defined earlier as a class of optically variable stars with emis- 
sion lines in the vicinity of bright or dark nebulosity by |Joy| 
( |1945| l. Ever since, the name-giving object of the class, TTau, 
has been regarded as the proto-type of young low-mass stars. It 
is a KO star with a mass of approximately 2 M© and an age of 
1-2 Myr, that has a modest and time-variable line of sight ex- 
tinction of Ay^lma g ( ,Kuhi. 1974,, Mel' nikov & Grankin|2005 



[Loinard et al.|2007l l. It resides in the Taurus molecular cloud at 
a distance of 148 pc ( |Loinard et al.|2007| ). 



* Based on observations obtained at ESO's La Silla/Paranal observa- 
tory under programmes 280.C-5035, 068.D-0537, 070.C-0544, 070.C- 
0162, 072.C-0593, 074.C-0699, 074.C-0396, 078.C-0386, 380.C-0179, 
382.C-0324, and 084.C-0783, as well as observations obtained at 
UKIRT. 



An "infrared companion", discovered by |Dyck et al.| ( [T982) l, 
is located approximately UH south of T Tau. It was soon realized 
that the bolometric luminosity of this source (T Tau S) rivals that 
of the optically visible star (T Tau N), and that it must be a self- 
luminous object rather than a dust cloud heated by the north- 
ern component. Even though TTauS is at least 10 mag fainter 



than the northern component in the optical (Stapelfeldt et al 



19981, it dominates the system flux at infrared wavelengths be- 
yond ^3 fim (Ghez et al.||1991| |Herbst et al.||1997t Beck et aL] 
2004), but it contributes only a minor fraction to the total flux at 



millimeter wavelengths ( ,Hogerheijde et al.,,1997; |Akeson et al. 
1998[ ). The southern component was itself found to be a close 



binary with components Sa and Sb at a projected separation of 
xiSOmas at the time of discovery (Koresko 2000), corresponding 
to only 7-8 AU at the system distance. 

T Tau S was shown to be strongly variable at near- and mid 
infrared wavelengths by Ghez et aL] (|1991)l, who detected a 
x2 mag increase in brightness at 2.2 /zm within ~ 1 yr, and a sim- 
ilarly large increase at 10 fim between two observations sepa- 
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rated by about five years. They attributed the brightness fluc- 
tuations to variable accretion and derived an accretion rate of 
3.6xlO-<'MQyr-i. 



Beck et al. ( 2004 1 presented a large set of (near-) simulta- 



neous K and L' photometry obtained over the course of about 
seven years starting 1995, and found TTauS to vary with a to- 
tal amplitude of nearly 3 mag at K-band and more than 2 mag at 
L' during this period. The fastest brightness variation they de- 
tected occurred in December 1997, when TTauS brightened by 
0.9 mag in the K band over only a seven day period. They ob- 
served the system during three consecutive nights in September 
1998 and another three nights in November 1999 in order to de- 
tect even faster variations, but did not detect any. The brightness 
variations in T Tau S were shown to be color-dependent in a K- 
L' vs. K color-magnitude diagram, in a "redder when fainter" 
fashion that closely resembles the ISM extinction curve (Rieke^ 
|& Lebofsky||1985! ). This behavior was not expected if the ob- 
served brightness fluctuations were due to variable accretion, 
in this case models predicted the opposite behavior, i.e. "bluer 
when fainter" (Ca lvet et al.|1 997 ), mostly because of the chang- 
ing contrast ratio between the blue stellar spectrum and the much 



redder disk spectrum as the accretion rate varies. Therefore Beck 
|et al.| ( [2004| l concluded that variable obscuration is the main re- 
sponsible factor for the observed large brightness variations in 
TTauS, although these authors did also find evidence for vari- 
able accretion in TTau S ( ,Beck et al.|2001[ l. 

A decade of spatially resolved observations of the southern 
binary pair have now covered a substantial fraction of the Sa-Sb 
orbit and have allowed direct determination of the orbital pa- 
rameters (|Kor esko 20 001 |Furlan et al.||2003t [Beck et al.HIOMl 
[Duchene et al.|2006t|Kohler et al.|2008 1. Current estimates still 
allow a substantial range for some of the orbital parameters, of 
which the orbital period and the mass of Sb are the most relevant 
for this paper Continued observations over the coming years are 
expected to yield a converged solution. The mass of the primary 
Sa is a!2.2 M0, the most recent estimates for the mass of Sb range 
from 0.4 to 0.8 Mq, and for the orbital period from a; 25 to ^lOO 
years (Kohler et al. 2008; Kohler 2008). The spatially resolved 



observations have also revealed Sb to show only minor bright- 
ness variations at 2.2 /zm (7?M5 *0.2 mag), while Sa has varied 
by more than 3 mag during the last decade at this wavelength. 

Here we revisit the physical mechanism causing the ob- 
served photometric variations of TTauS. This work was moti- 
vated by short-term mid-infrared variability that was serendipi- 
tously detected during an observing campaign with VLT/VISIR 
in February 2008, primarily aimed at characterizing the [Ne n] 
emission in the T Tau system. We reported on the [Ne 11] obser- 
vations in a separate paper ( van Boekel et al.|2009 1 and here fo- 
cus on the x26% increase in continuum brightness of TTau S at 
12.8 yum that we detected in two observations separated by only 
four days. In this paper we will distinguish between the very 
fast variations with modest amplitude as detected in our VISIR 
observations, and the long-term variations of much larger ampli- 
tude that are seen in the entire body of photometry taken over 
the past decades. These are not necessarily caused by the same 
physical mechanism. 

The paper is organized as follows. In Sect.|2]we describe the 
VISIR observations as well as additional previously unpublished 
infrared measurements performed with a range of facilities. The 
results of the infrared photometry are reported in Sect. |3] in 
which we also construct updated near- and mid-infrared light 
curves of the various components of the system. In Sect. |4. l| we 
explain why the short-term variability revealed in the VISIR ob- 



servations cannot be due to variable extinction and must instead 
be due to intrinsic luminosity variations, i.e. variable accretion. 
In Sect. 4.2 we address the nature of the long-term variability, 
which has previously been attributed to time-variable extinction. 
We present radiative transfer calculations of the disk of Sa that 
include variable accretion, and show that these models qualita- 
tively reproduce the observed long-term color-dependent near- 
and mid infrared brightness variations. In Sect. 4.3 we outline a 
tentative scenario for T Tau S in which the periods of high and 
variable brightness in which Sa has prevailed from a: 1990 to 
a2003 were due to enhanced accretion activity, quantitatively 
intermediate between an EXOR and a FUOR outburst, induced 
by gravitational perturbation of the Sa disk during the periastron 
passage of Sb in a 1995. We summarize our results in Sect. [5] 

2. Observations and data reduction 

The primary data set that motivated this paper comprises two 
12.8 /im images of the TTau system taken with VISIR at the 
VLT in February 2008. In addition, we have gathered an exten- 
sive set of photometric measurements from the literature, and 
present previously unpublished mid-infrared data obtained with 
the MAX instrument at UKIRT and the TIMMI2 instrument at 
the ESO 3.6m. We also present photometry of all K-band obser- 
vations performed with NACO at the VLT since 2001, including 
a number of previously unpublished data sets. 



2.1. VLT VISIR imaging at 12.8 /urn 



T Tau was observed with the mid infrared imager and spectro- 
graph VISIR dLagage et al. 2004| l, mounted on Melipal, the third 
of VLT's four 8.2 m Unit Telescopes. Images were taken through 
an approximately 0.23 yum-wide filter centered on 12.81 /i/m. 
The filter is centered on the [Ne II] 12.81 /vm line, but the emis- 
sion is dominated by the dust continuum and the neon line con- 
tributes only ^5% to the total flux integrated over the filter band- 
width. Observations were performed at two epochs, during the 
nights starting on 2008 February 3 and 7 (see Table[T]i. 

Standard chopping and nodding techniques were applied to 
deal with the high atmospheric and instrumental background in- 
herent to ground-based thermal IR observations. In a plain stack 
of the data the image quality of the first epoch is substantially 
worse than that of the second epoch. Closer inspection of the 
data cub^ shows that the position of the sources varied by sev- 
eral pixels between different chop half-cycles (here referred to 
as "frames"). In a sub-set of the individual frames the beams 
are strongly distorted, indicating that either the telescope or the 
chopping secondary mirror has moved during the integration. 
The remaining frames have the same beam quality as the data 
of the second epoch, where this problem was not encountered. 
Only frames with sharp, round beams were selected and com- 
bined after appropriate alignment. We combined the 25% best 
frames into our final images used for the analysis, and verified 
that our results are independent of this fraction as long as it is be- 
low a!60%. This reflects the image statistics of the first epoch, in 
which roughly 40% of the frames showed distorted beams. This 
gives us confidence that our method does not have unintended 
side effects on the photometry performed on the final image as 
long as all bad frames are removed. We applied the same proce- 
dure to the data of the first and second epochs. 



' The TTau system is very bright, allowing the beam shape in each 
individual frame to be assessed. Both T Tau N and T Tau S were detected 
at > 100 cr in every frame. 
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In order to determine the system sensitivity and the atmo- 
spheric transparency, photometric standards were observed im- 
mediately following the TTau observations during both nights. 
To correct for the airmass difference between the science obser- 
vations and the calibration measurements, synthetic atmospheric 
transmission profiles calculated with ATRAN ( |Lord|1992 l were 
integrated over the filter's spectral response and appropriate cor- 
rection factors were determined. Because the used narrow band 
Nell filter lies in a relatively "clean" region of the atmospheric 
transmission curve, the resulting airmass corrections are small 
(<6% during the first night, xl% during the second night). Any 
uncertainty in these corrections affects the absolute flux cali- 
bration but not the relative photometry between TTauN and 
TTauS. 

Fluxes of the individual components were obtained by per- 
forming PSF photometry, using the calibrators observed imme- 
diately after the science measurements as PSFs. We slightly 
smoothed the calibrator observed during the first night by con- 
volving it with a Gaussian of FWHM-Of.'l to correct for the 
relatively large difference in airmass between the science and 
calibration observation. 

The accuracy of our absolute photometric calibration is dom- 
inated by systematic effects, and is difficult to assess because 
we observed only one calibrator during each night. Thus, the 
stability of atmospheric transmission was not monitored, and to 
correct for diff'erences in airmass between the science and cali- 
bration observations we had to revert to theoretical calculations, 
as mentioned earlier. We adopted an uncertainty of 10% in the 
absolute calibration of our photometry. We emphasize that the 
relative photometry between T Tau N and T Tau S, i.e. of the N/S 
flux ratio, can be derived with a much higher precision of <2%, 
independent of absolute calibration uncertainties. 

2.2. Additional infrared observations 

2.2.1. UKIRT/MAX imaging 

The TTau system was imaged through a number of filters in 
the wavelength range 4.7-20 //m, on the nights starting 1995 
November 13, 1996 January 14 and 1996 August 27. The data 
of the latter epoch were presented earlier by Herbst et al. ( 1997 1. 
Here we will include the data in the 12.4 fim filter, closest to our 
continuum sampling wavelength of 12.8 fim, of all three epochs. 

Standard chopping and nodding techniques were applied and 
the reference star a Tau, used both for flux calibration and as 
PSF reference, was in all cases observed directly before or af- 
ter TTau, with similar sky conditions and airmasses. Separate 
frames from individual chop cycles (typical exposure times of 
0.2 s) were saved, and later combined using shift-and-add pro- 
cedures to obtain optimal image quality. Fluxes of the individ- 
ual components T Tau N and T Tau S were extracted with PSF 
photometry. The components of the southern binary Sa-Sb are 
not spatially separated in these observations. For further details 
with regard to the observing and analysis procedures we refer to 
IHerbstet al.| ( fT997] l. 

2.2.2. TIMMI2 imaging and spectroscopy 

The TTau system was observed with the TIMMI2 instrument 
( Reimann et al.|1998) mounted at the ESO 3.6 m telescope at La 
Silla observatory, Chile, on the nights starting 2002 February 2 
and December 24. Longslit spectra with a North-South slit ori- 
entation were taken, as well as imaging observations in the N8.9, 
N9.8, and N11.9 filters. Photometric and spectroscopic standard 



stars were observed regularly throughout the respective nights. 
The spectra have been previously published by Przygodda et al. 
(|2003)l, who did not attempt to extract the fluxes from TTauN 
and TTauS separately. In this paper, we have re-analyzed the 
whole data set, focussing on extracting the northern and south- 
ern component separately. 

The TIMMI2 images marginally resolve the N-S system 
{A/Dx(y:73 at 12.8 yum, compared to the a;0'.'69 N-S separation, 
and a 0'.'2 pixel scale). Individual fluxes for TTau N and TTau S 
were extracted from the imaging observations with PSF photom- 
etry. In order to correct for small differences in the actual point 
spread function between calibration and science observations, 
either the science or PSF reference observation was convolved 
with a Gaussian of FWHM<(y.'3, where the exact value was 
chosen such that the residuals in the fit were minimized. We kept 
as many parameters fixed as possible. In particular, we used the 
known positions of Sa and Sb with respect to N, and allowed 
only for a variable contribution of Sa and Sb, giving some free- 
dom to the photocenter of TTauS as a whole along the Sa-Sb 
separation. Thus, we have five free parameters in the fit: the (x,y) 
position of TTauN, multiplicative factors for N and S, and the 
relative contribution of Sa and Sb to the total flux of T Tau S. 

A similar procedure was applied to the longslit spectra, 
somewhat complicated because the 0'.'45 pixel scale in spectro- 
scopic mode undersamples the PSF. We modeled the observed 
two-dimensional spectra as recorded on the TIMMI2 detector 
with the minimum number of free parameters. We used an it- 
erative two-stage procedure aimed at minimizing the residuals 
between the observed and modeled signal. At each wavelength, 
the profile in the spatial direction was taken to be the sum of 
two delta functions with a fixed separation of 0'.'69, convolved 
with the profile of a calibrator star (stage 1). The free parame- 
ters are the amplitudes of both components, and the fit was done 
for all wavelengths sequentially. Then, in stage 2, the position 
and tilt of T Tau N on the TIMMI2 chip was varied, after which 
stage 1 was repeated, etc., until convergence was reached and 
the residuals were minimized. In short, we fitted the amplitudes 
of T Tau N and T Tau S at each wavelength separately, while the 
position of T Tau N was fitted to all wavelengths simultaneously 
and the N-S separation was kept fixed. 

2.2.3. VLT/NACO imaging at 2.2 //m 

Adaptive-optics-assisted imaging of the TTau system was per- 
formed with NACO, mounted on YEPUN, the fourth of VLT's 
8.2m Unit Telescopes, during numerous epochs since late 2001. 

Standard near-infrared data reduction methods were applied 
to the NACO images. They were sky-subtracted with a median 
sky image, divided by a flat field, and bad pixels were replaced 
by the median of the closest good neighbors. Finally, the images 
were visually inspected for any artifacts or residuals. 

We re-analyzed the entire data set, which includes five pre- 
viously unpublished observations. Relative photometry of all 
three stars was performed with T Tau N as the reference. A good 
model for the PSF is needed, which includes the seeing halo that 
is due to the imperfect adaptive optics correction. This requires 
some care, because the southern binary generally lies well within 
the seeing halo of the much brighter northern component. We 
constructed the PSF as follows. We first made an azimuthally av- 
eraged radial intensity profile I(R) by performing aperture pho- 
tometry with 128 circular apertures with logarithmically spaced 
radii, centered on TTauN, but with a 60° sector centered on 
T Tau S masked out (the resulting intensity profile is multiplied 
by 6/5 to correct for the masked out part). We found that at large 
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Fig. 1. Two-epoch VISIR imaging of the TTau system in a nar- 
row band filter centered on 12.81 jum. The left panel shows the 
system as it appeared on 2008 February 4, the middle panel as it 
did four days later, the right panel shows the difference. The fast 
brightening of T Tau S is apparent. All images are displayed in 
a linear stretch chosen in a way that TTauN peaks at the same 
level in both epochs. Values in the difference image have been 
multiplied by a factor of 4 for better display. Coordinates are in 
arcseconds relative to T Tau N. See section[3]for more details. 

radii, where only the seeing halo contributes, the intensity profile 
can be very well approximated with a quadratic fit in log(/(/?)) 
vs. R space. At radii of more than 22 pixels (a!0'.'292) we approx- 
imated the PSF with this fit. At smaller radii we used the actual 
image of T Tau N as our PSF. In this way we have a PSF that 
includes the inevitable AO artifacts close to the center and has 
a well approximated seeing halo, which is not contaminated by 
the southern binary. 

Higher-order AO residuals are seen in all NACO images of 
TTau and other bright sources. The most pronounced of these 
are four roughly static features at =sO'.'51 from TTauN located 
at »45, ~135, ~225, and »315 degrees East of North. The am- 
plitude of these features is approximately 1% of the peak flux of 
Sb, and they are not spatially coincident with TTauS. Thus our 
photometry of the Sa-Sb system is not compromised by our ap- 
proximation of the PSF with an azimuthally averaged fit at large 
radii, because any departures from the true PSF are negligible 
compared to the signal of the southern binary. 

The actual PSF photometry was performed with the 



Starfinder program (Diolaiti et al. 2000j), which finds point 
sources and determines their positions and flux multiplication 
factors in an image, given the PSF. We carefully checked the 
residuals, i.e. data minus the synthetic images provided by 
Starfinder, and found that at the position of Sb we get clean 
residuals, with occasional minor artifacts of <1% of the Sb peak 
flux. This shows that Sb is a perfect point source at 2.2 //m, at 
the resolution of NACO. Around Sa, however, we see substantial 
extended emission. This emission has the tendency to slightly 
boost the Sa flux in the PSF fitting because the PSF wings try to 
match the extended component. Consequently, the central peak 
emission is somewhat overestimated, resulting in slightly nega- 
tive residuals at the position of the Sa point source. We corrected 
for this by down-scaling the Sa fluxes, so that the residuals at 
the fitted position of Sa are non-negative. This required down- 
scaling the fluxes of Sa by typically 5% and always less than 
10%. See Sect. 3.2.3 and Fig.[3]for further discussion. 



3. Results 

3.1. Fast photometric variations of T Tau S at 12.8fim 

Figure[T] shows our 12.8 fim VISIR images of the TTau system 
taken during two epochs in early February 2008, separated by 
3.94 days. Though TTau shows relatively strong [Neii] emis- 
sion at 12.81 fim ( |van den Ancker et al.|1999| , the emission line 



contributes <5% to the total system flux integrated over the filter 
width, and our images are dominated by continuum dust emis- 
sion. 

The northern and southern components are spatially sepa- 
rated in both images and their fluxes could be determined indi- 
vidually. We measured the flux of TTauN to be 7.8 and 8.5 Jy 
during the first and second night, respectively. We attribute this 
difference to the uncertainty in our absolute calibration (see also 
then next paragraph), and in the further discussion we will as- 
sume the flux of T Tau N to be constant at our average value 
of 8.2 Jy, consistent with the flux of 8.3 Jy measured by |Ratzka| 
et al. '(2009') in December 2004. The fluxes of T Tau S were cor- 
respondingly scaled to 12.8 and 16.1 Jy during the first and sec- 
ond epoch, respectively (note that the values listed in Table[T]are 
not scaled but rather correspond to the directly measured fluxes, 
including absolute calibration uncertainties). 

3.1 .1 . How significant is tine detected briglitness cliange? 

Taking our measurements at face value, we observed a -1-9% 
(silcr) flux increase in TTauN and a -1-38% (*4cr) brightness 
increase in TTauS over the course of four days. The former 
is clearly not significant, whereas the latter is. The accuracy 
of these numbers is dominated by uncertainties in the abso- 
lute photometric calibration. The relative photometry between 
T Tau N and T Tau S is much more accurate, we conservatively 
estimate that we can determine this quantity to an accuracy of 
2%. Thus we can state, as a direct observational result, that the 
N/S brightness ratio changed by ^13cr (h-26%) in four days. 
Because T Tau N shows no significant brightness change in our 
observations, and it is known to be of approximately constant 
flux in the mid-IR, we have assumed that it did not change in 
brightness between 2008 February 4 and February 8. Thus, the 
flux increase in T Tau S becomes equal to the increase in flux ra- 
tio, i.e. +26% and a!l3cr. We adopt these values in the current 
analysis, and note that the assumption of a constant flux level for 
T Tau N has no qualitative and only a minor quantitative influ- 
ence on the reasoning and conclusions presented in this paper. 

3.2. Previously unpubiished infrared photometry 
3.2.1. UKIRT/IVIAX imaging 

The measured fluxes of TTauN and TTauS at 12.4 /im, ob- 
tained with UKIRT/MAX, are listed in Table |2] Within the un- 
certainties of the absolute flux calibration, T Tau N remained at 
a constant brightness whereas T Tau S varied substantially be- 
tween the three epochs of MAX observations. 



Table 1. Summary of the VLT/VISIR observations. All images 
were taken through the [Nell] filter. The last column lists the 
source brightness at 12.8 ;um. For the TTau observations, fluxes 
for the north and south components are individually given. The 
adopted fluxes of flie calibrators HD41047 and HD 75691 ai-e 



taken from 'C ohen et al.| ( |1999p . The absolute flux calibration is 
good to a; 10%. 



target 


observing epoch 


airmass 


^int 


F„(12.8 yum) 




(Universal Time) 




[s] 


N/S[Jy] 


TTau 


2008 Feb 04 02:09 


1.73 


143.5 


7.8/12.1 


HD 41047 


2008 Feb 04 02:42 


1.04 


54.3 


5.60 


TTau 


2008 Feb 08 00:40 


1.45 


137.8 


8.5/16.7 


HD 75691 


2008 Feb 08 01:17 


1.33 


52.0 


9.72 
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Table 2. UKIRT/MAX photometry of the TTau system at 
12.4 yum. 



date 



F,,(N) F„(S) 



1995 Nov 13 

1996 Jan 14 
1996 Aug 27 



9.4 
9.8 
9.7 



16.3 
19.7 
15.8 



3.2.2. TIMMI2 imaging and spectroscopy 

Figure[2]shows the results of our PSF photometry and PSF spec- 
troscopy performed on the 2002 TIMMI2 data. The spectra were 
rebinned to a resolution of s»25, and the error bars indicate the 
standard deviation within each bin, i.e. they reflect statistical 
fluctuations but not necessarily systematic uncertainties in, e.g., 
the telluric calibration or errors arising in separating the fluxes 
of both components in the marginally resolved and undersam- 
pled data. The "gap" in the December spectra around 9.2 yum is 
due to a defect detector channel. The spectra taken in February 
and December 2002 were scaled by factors of 1.05 and 1.29, 
respectively, to match the flux levels derived from the imaging 
observations. Such factors may arise from uncertainties in the 
absolute calibration of the spectra, i.e. due to diflFerent slit losses 
between calibration and science observations. The photometry 
was not scaled. 

The spectra of TTauN taken in February and December 
2002 are identical in shape within uncertainties and show a weak 
silicate emi s sion feature that has been reported earlier (e.g. Ghez 
T99T] [Herbst et al.|[T997l [Ratzka et al.|[2009] ). In partic 



et al. 



ular the February spectrum of T Tau S clearly shows the well 
known silicate absorption feature. The February spectrum shows 
TTauS to be approximately as bright as TTauN in the con- 
tinuum redward of the silicate feature. This spectrum is nearly 
identical to that measured by Ratzka et al.| ( 2009] l. In December 
the shape of the TTauS spectrum has changed substantially: 
whereas at wavelengths below ^10 jum it is still very similar to 
the February spectrum, beyond 10 yum the red wing of the sili- 
cate absorption feature no longer rises stongly, but instead has 
flattened considerably. 

We remind the reader that in our TIMMI2 observations 
T Tau N and T Tau S are only marginally spatially separated, due 
to the the modest primary mirror diameter of the ESO 3.6m tele- 
scope. This limits the accuracy with which the fluxes of the in- 
dividual components can be extracted. However, with our robust 
approach using a minimum of free parameters (see Sect. 2.2.2| i 
we get a fair agreement between the imaging and spectroscopic 
observations, which is reassuring. 

3.2.3. VLT/NACO imaging at 2.2 jum 

In Table[3]we summarize the VLT/NACO K-band photometry. 
Most of these data were previously published, but with an em- 
phasis on the relative astrometry of the various components 
rather than the photometry (e.g. |Kohler et al.||2008] [Brandner] 
et al. 20021. Here we present the photometry of all existing 
NACO measurements, including five previously unpublished ob- 
servations. We did not attempt to perform an absolute photomet- 
ric calibration. Instead, all photometry of Sa and Sb is relative 
to TTauN. The formal errors on the photometry as given by 
Starfinder are generally less than 1% for Sb and less than 2% 
for Sa. However, given the non-negligible residuals, in particular 
surrounding Sa, these estimates seemed too optimistic. We as- 
signed errors of 0.02 mag to the photometry on Sb and 0.05 mag 
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Fig. 2. PSF photometry (squares) and PSF spectroscopy per- 
formed with TIMMI2 in 2002. The northern and southern com- 
ponents are spatially marginally resolved (» Q.9A/D at 12.8 fim). 
See Sect. l2.2.2l for details. 



on that of Sa. While these estimates may be somewhat conser- 
vative, the amplitude of the photometric variations that we are 
discussing here is vastly larger than the uncertainties, and thus 
our simple error estimates serve the current cause well. 



As described in Sect. 2.2.3 we found substantial spatially ex- 
tended emission around Sa. Around Sb, no such emission was 
found. In Fig.|3]we show the average residuals of all observa- 
tions, aligned in a way that Sa is always at the center of the 
image as indicated with a + sign. The positions of Sb, one for 
each epoch, are also indicated with + signs. The residual emis- 
sion surrounds Sa, but shows a central cavity. Sa is not in the 
center of this cavity, but rather in the north-west corner The 
emission level of the brightest parts of the extended emission 
corresponds to x5% of the peak flux of Sb. The emission has its 
largest extent towards the south. Note that Sa itself has moved 
by approximately 24 mas (assuming a mass ratio of 2.2:0.6 for 
Sa:Sb), or about 1/3 of the diameter of the "central cavity" dur- 
ing the period covered by our NACO observations. The spatially 
extended material we see need not necessarily move along with 
Sa. 

One may ask whether the cavity seen is real or an artefact of 
imperfect subtraction of the Sa point source. After careful con- 
sideration, we conclude that it must be real for the following 
reaons: (1) the spatially extended emission and the central cav- 
ity are persistent over time; (2) the integration times were chosen 
so that T Tau N, which is used as PSF reference, stays below the 
linearity limit of the NACO detector; (3) the cavity is offset from 
the Sa point source (see the dashed contours in Fig.|3]i; and (4) 
during the observations of the last epoch (October 2009) the field 
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Fig. 3. Spatially extended K-band emission around Sa. Shown is 
a contour plot of the NACO images after subtraction of the best- 
fit point sources for TTauN, Sa, and Sb (see Sections [2.2.3 and 
3.2.3 for details). All epochs were aligned and centered on Sa. -H 
signs mark the positions of Sa and Sb. The inset shows the PSF 
with the 13.26 mas pixels of the sl3 camera on the same scale as 
the main plot. The dashed contours show the best-fit point source 
of Sa, drawn at 25% and 75% of the Sa peak flux. 



was rotated by a; 130°. This leaves the orientation of PSF on the 
detector unchanged, but rotates the sky. The extended emission 
and the central cavity are detected at the same location during 
this epoch, even though these regions of the sky now correspond 
to a completely different part of the PSF. This proves that the 
central cavity cannot be an artifact due to a subtle asymmetry in 
the PSF 

Our aim here is not to assess the nature of the extended emis- 
sion in detail. The main point we wish to stress is that there is 
material around Sa on scales of ~5 to 20-30 AU. This scale is 
larger than the radius of the Sa disk, which must be <5 AU, but 
substantially smaller than the extended H2 2.12 /vm line emis- 
sion seen in IFU spectra by |Beck et al.| (j2008 ) and |Gustafsson| 
[etaL] ( |2008] l on scales of ^50-200 AU. The H2 IFU data do not 
properly resolve the scales on which we see the extended emi- 
sison in the NACO data, but they suggest the H2 emission to be 
very faint close to Sa. This argues for K-band continuum radi- 
ation from Sa scattered by dust as the most likely explanation 
for the emission we see. Whether the "central cavity" is truly 
devoid of material, or whether there is material present at these 
locations which has no direct line-of-sight to Sa and the near-IR 
bright inner regions of its disk and hence does not show up in 
scattered light, remains to be investigated. 

3.3. Which components vary by how much? 

The photometric observations of the TTau system performed 
over the last decades, which form the basis of the current and 
previous variability studies, had a range of spatial resolutions 
depending on the facility used and the wavelength of observa- 
tion. At near-infrared wavelengths, all observations used here 
spatially resolve T Tau S from T Tau N, but only the more recent 



Table 3. Summary of the VLT/NACO photometry in the Ks 
band obtained between 2001 December and 2009 October All 
magnitudes are relative to T Tau N. 



date 

2001 Dec r 

2002 Dec 15 

2003 Dec 12 

2004 Dec 9 

2006 Oct 11 

2007 Sep 16 

2008 Feb 1 
2008 Oct 17 

2008 Nov 6 

2009 Feb 18 
2009 Oct 7 



AA:^ (Sa-N) 
1.67±0.05 
3.61±0.05 
4.75±0.05 
4.79+0.05 
3.14+0.05 
3.09+0.05 
3.07+0.05 
3.46+0.05 
3.84+0.05 
3.45+0.05 
3.05+0.05 



AKs (Sb-N) 
3.02+0.02 
3.14+0.02 
3.25+0.02 
3.31+0.02 
2.94+0.02 
2.92+0.02 
3.05+0.02 
3.07+0.02 
2.88+0.02 
2.88+0.02 
3.07+0.02 



Table 4. References for infrared photometry used in Figs.|4j[5j 
and [6] The first three columns give the reference to the origi- 
nal work, the wavelength of observation, and the facility used, 
respectively. In the last two columns we indicate whether the 
observations spatially separated the northern and southern com- 
ponent and the Sa-Sb sub-components of TTauS, respectively. 



original work 



facility 



N-S Sa-Sb 
spat, resolved? 



Beck etal. (2004) 
Duchene et al. ( 2005| 
Herbst et al. (20071+ 
iGhez et al.|(|l991|) 
Ratzka et al (2009) 
v.d.Ancker etal. (1999) 
this work 

Ratzka et al. (2009) 
this work 
this work 
this work 



+ 



+ 



K, L' 
K, L' 
Ks,L' 
10.1, 12.5 
12.8 
12.8 
Ks 
12.8 
12.4 
12.8 

12^8 

^Includes also data from WIRO, Keck, and Gemini North. The Sa-Sb pair 
remained unresolved in most of these observations, ut several spatially resolved 
data points are included as well. Besides a large number of their own measure- 
ments. Beck et al. include literature data from Dyck et al.' ('1982X 'Beckwith] 
et al. 1 19841, Maihara & Kataza 1 1991 1, Ghez et al. ( 1991 1, Tessier et al. ( 1993|7 
Kobayashi et al. ( 19941, Simon et al. 1 19961, Roddier et al. (20001, Herbst et al. 



IRTFt 
Keck 
VLT 
Hale 5m 
VLTI 
ISO 
VLT 
Spitzer 
UKIRT 
ESO 3.6m 
VLT 



-I- 



_t 



( 19971, White & Ghez (2001 1, Koresko (20001, Kasper et al. (20021, Duchene 
et al.|^2002^ , and |Furlan et al.,^2003j . 

"''We applied scaling factors of AK=-0.18mag and AL'=+0.32mag to all 
components, so that the magnitudes of T Tau N match the mean literature values 
of K=5.52and L'=4.32. 



observations resolve the southern Sa-Sb pair At mid-infrared 
wavelengths, most observations spatially resolve TTauS from 
T Tau N, the space-based measurements yield only the cumula- 
tive flux of the whole system, and some measurements spatially 
resolve the Sa-Sb binary using special techniques. 

If one component of the triplet is the dominant source of vari- 
ability, we can with some care use all measurements to constiTict 
a light curve of this source, including those observations that do 
not spatially separate all three components. We will argue that 
TTauSa is the main responsible for the photometric variations 
of the whole system at mid-infrared wavelengths, and that it is 
dominates the variability of T Tau S in the near-infrared, as was 
previously done by e.g. Duchene et al. (2005i 

TTauN. Historic optical photographic plates have revealed 
that the optically visible component T Tau N was varying irreg- 
ularly, rapidly, and strongly between 1858 (beginning of data 
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taking) and -1917 ( |Lozinskii| 1949] [Beck & Simon|200l| . Since 
then, the optical brigthness of T Tau N has shown comparatively 
minor (<1 mag) deviations from its average value of B^ll mag, 
which is approximately the maximum brightness reached during 
the period of irregular variations prior to 1917, except for brief 
periods in 1925 and 1931 when the star was substantially fainter 
These variations have been attibuted to time-variable line-of- 
sight extinction caused by dynamic, structured dust clouds pass- 
ing in front of T Tau N, possibly related to gravitational inter- 
action between the stars and circumstellar matter in the system, 
or out-flowing material from T Tau S ( [Herbst et al 



1994 



Beck 



|& S imon 2001). Such variations in foreground extinction may 
cause infrared variability as well, but the amplitude is expected 
to be much smaller, in particular in the mid-infrared. 

Since the beginning of near- and mid-infrared observations 
in the 1980s, TTauN has not shown substantial variations at 



these wavelengths (for the near-lR see Beck et al. 20041. At 



12.8 /vm, the fluxes reported in the literature and the new data 
presented here average around 8.1 Jy, with a standard deviation 
of 0.9 Jy ( |Ghez et al.||l991| jSkemer et al.| |2008; Ratzk a et al. 
2009 1. Considering the limited accuracy of ground-based ab- 
solute photometry in the mid-lR and that these measurements 
were taken through different filters, requiring some extrapola- 
tion to our sampling wavelength of 12.8 /im, we conclude that 
the existing measurements show no evidence for significant mid- 
IR variability of T Tau N. Minor variations with an amplitude of 
<10% cannot be excluded based on the currently available data. 

TTauS. Near-infrared imaging on AO-assisted 10m class 
telescopes performed since the early 2000s has revealed Sa to 
be variable by over 3 mag at 2.2 jum, whereas Sb shows only 
modest variations with an RMS of 0.2 mag around a mean value 
of Ks~8.6 (Note that Koresko (2000 1 finds a substantially lower 
flux of K=9.37±0.25 for Sb in speckle holographic imaging per- 
formed in late 1997, assuming T Tau N has K=5.52). At 12.8 fj.m 
T Tau S has varied in brightness from ^5 to ^27 Jy between the 
epochs covered ( jGhez et al.|1991 ). Recent N-band observations 
that use mid-IR adaptive optics and interferometric techniques 
have spatially resolved the Sa-Sb pair in the N-ban d, and showed 
that Sb has a brightness of <2.5Jy at 12.8 fim (Skemer et al. 
Ratzka et al.|2009| l. Thus, it appears very unlikely that Sb 



2008 



contributes dominantly to the total mid-infrared flux of T Tau S 
at any epoch, and thus also not to its variability. This is agrees 
completely with the absence of large near-infrared variations in 
Sb. As argued in the previous paragraph, T Tau N shows no sig- 
nificant IR variability. Therefore we can reasonably attribute the 
vast majority of the near-infrared variability of TTauS and the 
mid-infrared variability of the whole TTau system to Sa only. 

We can now deduce the flux of Sa from observations that do 
not spatially separate all components of the system. In the near- 
infrared, where all used measurements spatially resolve T Tau S 
from T Tau N, we can obtain the magnitudes of Sa alone by tak- 
ing Sb to be of constant brightness at K=8.6 mag ( [Duchene et al.| 
|2005, this work). For epochs during which TTauS was brighter 
than K=8.0mag, the 0.2 mag RMS variations of Sb will intro- 
duce an error of <0.3 mag on the thus derived brightness of 
Sa. Because we are interested in global trends rather than high- 
precision photometry of Sa, uncertainties of <0.3 mag are ad- 
equate for our purposes (the total variations we are studying 
amount to AK>3 mag). We conservatively disregard any near- 
infrared observations in which the Sa-Sb pair is not spatially sep- 
arated and the total brightness of T Tau S is less than K=8.0 mag 
in the [Beck et aL| ( |2004[ ) data. Likewise, we calculate the L' 
brightness of Sa from the total brightness of TTauS assuming 



that Sb has a constant brightness of U-6.25 mag ( Duchene et al. 
2005). Because Sa is substantially brighter than Sb in L', the 
uncertainties on the derived L' magnitudes of Sa will be smaller, 
though strictly speaking we do not know how variable Sb is in L' 
because no multi-epoch, spatially resolved observations in this 
band are available. 

In a similar fashion, we can infer the brightness of Sa at 
12.8 from observations that do not spatially separate Sa from 
Sb (as is the case for most ground-based observations used here), 
or from space-based observations that yield only the cumula- 
tive flux of the entire T Tau system. T Tau S was much brighter 
than Sb at all covered epochs, and T Tau S was approximately 
as bright as (Spitzer epoch) or substantially brighter than (ISO 
epoch) TTauN in the space-based observations. Therefore the 
uncertainties in the derived 12.8 yum flux for Sa are <0.5 mag 
for the epochs during which Sa was faintest and substantially 
smaller during those epochs when Sa was bright. 

In summary, we attribute the vast majority of the observed 
near-infrared photometric variations in TTauS and the mid- 
infrared variability of the entire TTau system to Sa only. We 
acknowledge the minor near-infrared variations detected in Sb 
and potential <10% mid-IR variability in TTauN, but these are 
negligible compared to the Sa variability and have no qualitative 
and only a very minor quantitative influence on our discussion. 

3.4. Updated light curves of the TTau system 

Figure|4] shows a compilation of photometric measurements of 
the TTau system between 2 and 13 jum, spanning a time inter- 
val of nearly three decades. The data shown consist mostly of 
available literature photometry, supplemented by the new mea- 
surements presented in this work. Similar curves showing a sub- 
set of these data have been presented earlier ( |Ghez et al.|1991[ 
Beck et al.|2004) ). The references for all data shown are given in 



Table|4| 

For the N-band fluxes (10.1 fim broad band fluxes for the 
earliest epochs, continuum fluxes at or near 12.8 fj.m for the 
epochs since the early 1990s) we show both the total system 
flux (N+Sa-HSb) and the flux of TTauS (Sa+Sb) only. The N- 
band fluxes were converted to magnitudes using zeropoints of 
41.0Jy at 10.1 fim and 25.6 Jy at 12.8 yum. Note that only the 
fluxes at 12.8 yum are used in the further analysis. For the L' 
magnitudes we separately show the total flux of T Tau S (Sa-i-Sb) 
and the brightness of Sa alone. To obtain the fluxes of Sa from 
measurements that did not spatially resolve the Sa-Sb pair (Beck 
|et al.|[2004| see Table[4]i we assumed that Sb was constant at 
L' =6. 25 mag. In two epochs of AO-assisted imaging the L' 
fluxes of Sa and Sb could be determined directly (Duchene et al. 
|2005HHerbst et al.|2007j ). Also for the K-band we show the total 
flux of TTauS (Sa-i-Sb) and the flux from Sa alone separately. 
Besides a large number of spatially unresolved measurements, 
more than a dozen AO-assisted observations are shown that spa- 
tially separate the Sa-Sb pair. 

3.5. The color-magnitude behavior of the long-term variability 

The multi-wavelength infrared light curve presented in Fig.|4]al- 
lows us to investigate the relation between infrared brightness 
and SED shape of the infrared companion. The first comprehen- 
sive study of this kind was presented by |Beck et al.| (|2004), who 
presented a K-L' vs. K color-magnitude diagram of T Tau S and 
showed that its near-infrared color varies with brightness in ac- 
cordance to the ISM extinction law. In Fig.|5]we show a CMD 



R. van Boekel et al.: Variable accretion as a mechanism for brightness variations in T Tau S 

T Tau system light curves 



-1 1 1 1 1 r 

red: north + south 
block: south (So + 
green: So 
orange: Sb 



-1 1 1 1 n 1 1 1 1 n — n — n r- 

12.8 /J-m'^ 4i" -i- -i- 



T T T- 




— 



2 — 



01 



10 



12 



1980 



1985 



2005 



2010 



Fig. 4. Updated light curves of various components of the TTau system at various infrared wavelengths. References for all mea- 
surements shown are given in Table|4] At the top (roughly between and 3 magnitudes) are N-band measurements. Two epochs 
in which the observations were space based and did not spatially resolve the N-S pair are indicated with '1' sub-scripts. The two 
epochs of TIMMI2 observations that only very marginally resolve the N-S pair are marked with '2' sub-scripts. In red we plot 
the total system brightness (N+S), in black we plot the magnitude of the TTauS (sum of Sa and Sb). In the middle of the graph 
(roughly between 3 and 6 magnitudes) we show L' magnitudes of TTauS. The black diamonds show the total flux of TTauS, the 
green diamonds show the flux of Sa only. The observations of the last two epochs (large diamonds) spatially resolved Sa and Sb 
( Duchene et al.|2005 Herbst et al.|2007 1, here the magnitude of Sb is indicated in orange. At the bottom of the graph (roughly at 6 
magnitudes and below) we show K-band magnitudes of T Tau S. Triangles indicate measurements in which the Sa-Sb pair was not 
spatially resolved, green and orange squares indicate the magnitudes of Sa and Sb, respectively, deduced from spatially resolved 
observations (Sa-i-Sb in black). 



similar to that of Beck et al., except that we plot the magnitudes 
of Sa only (the Beck et al. diagram shows TTauS, i.e. the sum 

of Sa and Sb). 

In Fig.[6]we extend the study of Beck et al.| ( |2004 ^ to longer 
wavelengths by compiling a K-[12.8] vs. K color-magnitude di- 
agram of Sa. The 2.2 /im and 12.8 /im data were generally not 
obtained simultaneously and have time lags of up to nearly two 
months. Hence this discussion is necessarily limited to the long- 
term variability. Because Sa has been observed to vary substan- 
tially on timescales of days, we assigned uncertainties to the data 
points to account for the short-term variability in an approxi- 
mate sense as follows. We took the fastest observed variations at 
12.8 yum (0.26 mag in four days, this work) and scaled this with 
the square root of the time lag between the 2.2 and the 12.8 /im 
observations. This choice of scaling is somewhat arbitrary but 
reflects the apparent randomness of the observed short-term fluc- 
tuations. 

4. Discussion 

In this section we will first discuss why the fast mid-infrared 
variability that we detected excludes variable extinction as a vi- 
able mechanism for the observed short-term mid-IR brightness 
fluctuations in TTauS, and argue that these must instead be due 
to variations in intrinsic luminosity. Then we will present a ra- 



diative transfer disk model of Sa, including both an active ac- 
cretion and a passive reprocessing component, that qualitatively 
reproduces the observed long-term color-magnitude variations 
when the accretion rate is varied. Thus we will show that vari- 
able accretion may also explain the long-term variability, but do 
not exclude variable extinction as the prime mechanism respon- 
sible for the long-term variations. Lastly, we sketch a tentative 
scenario for T Tau S, in which the » 15 yr period of irregular and 
enhanced brightness that we witnessed in the recent past is in- 
duced by the periastron passage of Sb, gravitationally disturbing 
the Sa disk and triggering an accretion outburst that may be rem- 
iniscent of EXOR variables. 



4.1. Can the observed short-term brightness fluctuations be 
caused by variable extinction? 

For variable extinction to be a viable cause of the observed 
brightness fluctuations, it is a necessary condition that the 
opaque medium causing the extinction can (un-) cover the emit- 
ting region within the timescale of the variations. Because only 
Sa shows large brightness variations whereas TTauN and Sb 
do not, one may assume the absorbing "screen" to be local to 
the system, i.e. at a distance of a!l48pc. Thus we can estimate 
the minimum speed at which the screen must travel to cause the 
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observed variations. To this purpose, we will first calculate the 
minimum size of the emitting region. 

Our sole constraint is that TTauS brightened by +3.3 Jy 
within 94.5 hours at 12.8 jim. Because at 12.8 yum the emission 
is completely dominated by thermal emission from dust grains, 
the maximum intensity the emitting region can have is that of a 
~1 500 K black body (at higher temperatures, the dust evaporates 
and the material loses its IR opacity). Thus the minimum solid 
angle of the region that needs to be (un-) covered is: 



Q = Fy/By(l500K) 



1.9x10"'^ sr. 



where Fy - 3.3 iy, and By denotes the Planck function (evalu- 
ated at 12.8 jum). At the system distance of 148 pc, and approx- 
imating the emitting region by a circl^ this yields a diameter 
of a;0.48 AU. The minimum required speed for the absorbing 
screen to (un-) cover this region within the available four days 
then becomes ^210kms~'. 

This velocity is much higher than the velocities one may ex- 
pect in the close environment of TTauSa. The Kepler speed 
at 0.24 AU, the absolute minimum distance from the star at 
which the screen may be positioned in order to cover a region of 
0.48 AU in diameter, is si90kms"\ assuming a mass of 2.2 M© 
for Sa ( [Kohler et al";][2008| [K5hier| [2008). This velocity corre- 
sponds to a circular orbit, screens orbiting on eccentric orbits 
may reach velocities of up to ^130kms ' at this distance. If 
the absorbing screen would be on an eccentric orbit, though, it 
would itself become warmer and brighter at 12.8 fim on closest 
approach, thus at least partially compensating for its own dim- 
ming effect. Thus, extinction caused by a "structure" existing 
within the disk, e.g. a warp or spiral arm, cannot be responsible 
for the observed IR variability. Moreover, the velocity disper- 
sion within a molecular cloud is only a few kms ' (e.g. iLarson 
1981| l, which excludes dust clouds in relatively close vicinity of 
the star, but not directly related to it, as viable absorbing screens. 

In the above analysis we have made several simplifications, 
which were all chosen to lower the required velocities for the 
observing screen, i.e. to favor the variable extinction scenario. 
The goal was to show that variable extinction can be ruled out, 
even with these unrealistically conservative assumptions: 

• The required velocity of the obscuring screen that we de- 
rived is the absolute minimum value that does not violate 
fundamental laws of physics (Planck's or Kepler's law), un- 
der the assumption that the observed IR continuum radiation 
is thermal dust emission. However, in reality, the region of 
the disk that emits the bulk of the flux at 12.8 //m is likely 
much larger, leading to a correspondingly higher minimum 
required velocity. For a typical disk model around an object 
with approximately the stellar parameters of TTauSa, the 
"size" of the disk (here defined as the region within which 
75% of the flux at 12.8 /urn is emitted) is about 15 AU in 
radius (e.g. van Boekel et al. 2005 1. In the special case of 
T Tau Sa, where the disk is likely truncated on the outside at 
a radius of 3-5 AU due to tidal interactions with TTauSb, 
the region we see at 12.8/im probably encompasses the en- 
tire disk. Conservatively assuming an outer radius of 3 AU 
for the disk, an absorbing screen would need to travel at 
«700kms"' in order to (un-) cover 26% of the emitting re- 
gion within four days, whereas the Kepler velocity at this 
radius is only a;26kms"'. 



• We assumed the absorbing screen to be 100% opaque at 
12.8 fj.m. If the screen's transparency at this wavelength is 
considerable, a larger area needs to be (un-) covered within 
the available time and the required velocity becomes corre- 
spondingly higher 

We conclude that variable extinction cannot cause the short-term 
brightness fluctuations observed at 12.8 yum. The absolute min- 
imum velocity required for an absorbing screen to (un-) cover 
the emitting region within the available time is about a factor 
of 2 higher than the highest possible speed such a screen could 
have at the relevant location. A more realistic estimate yields a 
discrepancy of a factor of a!20-30 between the required and the 
physically possible velocity. 

Note that if we apply the same simple model to the fastest 
variations detected by |Beck et al.| ( [2004) l at 2.2 fim (^0.9 Jy in 
seven days|^ we cannot exclude variable extinction: a 1500 K 
blackbody emitting the corresponding flux would have a di- 
ameter of ^0.15 AU and the minimum speed needed for an 
opaque screen to (un-) cover this region within 7 days is 
only a!37kms"'. This is substantially below the correspond- 
ing Kepler speed of a!l60kms and thus variable extinction 
is a perfectly viable explanation for the observed variations. 
Considering that the stellar photosphere may contribute signifi- 
cantly to the total flux at 2.2 jjm, the case for variable extinction 
becomes even stronger Only the long wavelength of our VISIR 
observations allows rejection of the variable extinction hypoth- 
esis. 



4.2. Can the observed long-term brightness fluctuations be 
caused by variable accretion? 

In the previous section we showed that apparent brightness fluc- 
tuations due to variable extinction cannot be the cause of the 
observed short-term mid-IR photometric variations of TTauSa. 
Therefore the variability must be due to intrinsic changes in the 
luminosity of Sa. A high and strongly variable accretion rate 
would naturally lead to such luminosity changes in a young star 



^ The brightness of TTauS increases from K=7.4mag on 1997 
December 6 to K=6.5 mag seven days later jBeck et al.| |2004| l. 
Assuming Sb was constant at K=8.6 mag, this implies that Sa went from 
K=7.84mag to K=6.67 mag, i.e. it brightened by a; 1.2 mag, or kO.9 Jy. 



Table 5. Summary of the parameters of the radiative transfer 
model described in Sect. |4.2| and Figs.|5]|6j and|7] 



^ The argument is qualitatively independent of the exact shape of the 
emitting region and we adopt the simplest possible geometry. 



parameter 


symbol 


value 


unit 


stellar radius 


R, 


3.5 


Ro 


stellar mass 


M, 


2.2 


Mo 


effective temperature 




5500 


K 


stellar luminosity 


L. 


10 


Lo 


disk mass 




10-3 


Mo 


exponent surface density 


p 


-1 




inner radius passive disk 


Rm 


0.5 


AU 


outer radius passive disk 


Rout 


5.0 


AU 


scale height at inner edge 


Ho 


0.021 


AU 


flaring index 


1 


2/7 




accretion rate 


Mac, 


lo-^ lo-^ 10-^ 
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Fig. 5. K-L' vs. K color-magnitude diagram similar to that pre- 
sented by Beck et aT1 ( |2004| l, except that we plot only Sa ( |Beck| 
et al. ( 2004[ l plot the sum of Sa and Sb). The diamonds repre- 



sent indirect measurements calculated from the data by Beck 
|et al.| ( [2004l l, in which Sa-Sb remained spatially unresolved, as 
suming that Sb was of constant brightness at K=8.6mag and 
L' =6. 25 mag, respectively. The diamonds in the upper left part 
of the diagram, in a darker shade, are less reliable because Sa 
was comparatively faint during the epochs of these measure- 
ments. The cross and triangle are from'Duchene et al. ( 2005 [ l and 
Herbst et al.| ( j2007 ), respectively, whose measurements did spa- 
tially resolve the Sa-Sb pair Over-plotted are radiative transfer 
model magnitudes, labe led with their respective accretion rates 
in Moyr"\ see Sect.lOand Tablels] 



Fig. 6. K-12.8 vs. K color-magnitude diagram of TTauSa. The 
diamonds represent measurements in which the Sa-Sb pair was 
spatially resolved in the K-band, crosses denote measurements 
in which we calculate the K-band magnitude of Sa from the to- 
tal flux of TTauS assuming Sb has K=8.6. The N-S pair was 
spatially resolved at 12.8 /im, except for the two points encir- 
cled in grey, for which the 12.8 fim flux was calculated from 
the total flux of T Tau assuming T Tau N has a flux of 8.2 Jy. The 
point surrounded by a grey box represents a broad band measure- 
ment centered at 12.5 /vm by Ghez et al. ( |1991 1. The uncertain- 
ties arise mostly because the near- and mid-infrared observations 
were not performed simultaneously (see Sect. 3.5 1. Over-plotted 
are radiative transfer model magnitudes, la bele d with their re- 
spective accretion rates in Moyr"', see Sect. 4.2 and Tablels] 



like Sa, as already argued by Ghez et al. (19911. Thus, vari- 
able accretion is the prime candidate mechanism to explain the 
observed short-term variability. One may ask whether it can 
explain the observed long-term photometric variations, which 
are of much larger amplitude, as well. The variable accretion 
mechanism has previously been deemed unfit to explain the 
"bluer when brighter" color-magnitude behavior of the long- 
term brightness fluctuations ( Beck etal. [200 4 ). In this section we 
will show that constructing a disk model that qualitatively repro- 
duces the observed color-magnitude behavior is in fact possible, 
with parameters that remain in a realistic range. 

Our model of the Sa disk consists of two components: an 
"active" inner disk in which luminosity is generated by the re- 
lease of gravitational energy from accreting material, and a "pas- 
sive" outer disk that is irradiated by the central star and the in- 
ner disk, and reprocesses the absorbed optical/NlR photons into 
longer wavelength infrared radiation. The model parameters are 
briefly discussed here and are also listed in Table |5] They are 
chosen so that they represent Sa and its disk to the best of cur- 
rent knowledge. In particular the small outer radius of the disk 
and the near edge-on orientation are key features (e.g. |Duchene 
let al.|2005j l. 

The active, inner part is modeled with an optically thick, ge- 
ometrically thin accretion disk with a radial temperature pro- 
file following the |Shakura & Sunyaev] ( |1973| l model. In addi- 
tion, a "hot-spot" covering 1% of the star, representing the ac- 
cretion shock of infalling material hitting the stellar surface, 
is included by adding a blackbody of the corresponding tem- 



perature and luminosity to the stellar spectrum. For the pas- 
sive, outer disk we used the 2D radiative transfer code radmc 
(DuUemond & Dominik 2004) to calculate the temperature 



structure and emergent infrared SED. The passive disk part has 
an inner radius of Ri^ - 0.5 AU, an outer radius of /?out = 5 AU, 
a mass of M^isk = 10"^ Mq (assuming a gas to dust ratio of 
100) with a radial surface density profile of Zoc/? a flared 
geometry with a flaring inde?([^ of y = 2/7 and a scale height 
of i/o = 0.021 AU at the inner disk edge of the passive, re- 
processing disk {HjR - 0.042, chosen to be the value corre- 
sponding to hydrostatic equilibrium in the "high" state). The 
outer radius is much smaller than that of typical circumstel- 
lar disks because of tidal truncation by Sb, whose highly ec- 
centric orbit brings it to within Ril2 AU from Sa at closest ap- 
proach. This naturally explains the non-detection of the Sa disk 
at mm wavelengths by | Hogerheijde et al.|(1997) and Akeson 
jetal]([T998|, as argued before by |Koresko| ( |2000| l 

By varying the accretion rate, we change the amount of en- 
ergy that is released in the center of the system, heating the cir- 
cumstellar material. The luminosity of the stellar photosphere 
dominates at low accretion rates, providing a "base level" for 
the energy output. At accretion rates above a!7xl0"^ Moyr"' ac- 
cretion becomes the main energy source and the heating rate of 
the circumstellar material scales roughly linearly with the accre- 
tion rate. Variations in the central luminosity lead to essentially 



* The flaring index y is defined so that the ratio of the disk scale 
height H and the distance to the central star R scales like H/R oc R^. 
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seen at an inclination of 74° and with a foreground extinction of 
15 mag agrees fairly well with the observations (the extinction 
matches the value of the foreground extinction observed toward 
Sb (Duc hene et al.|2005] )). This particular inclination yields the 
best fit for the set of model parameters that we chose, but other 
choices of disk parameters will result in slightly different best- 
fit inclinations. Thus, our model does not accurately constrain 
the inclination of the Sa disk, but we stress that a near edge-on 
inclination is required to obtain the correct color-magnitude be- 
havior of the variations, regardless of the specific choice of disk 
parameters. 

In Fig.|7]we show the model SEDs for accretion rates of 10 
10"^, and 10"^ Moyr"', as well as the observed infrared pho- 
tometry. We also show the model for the highest accretion rate 
as it would appear if the foreground extinction is omitted (grey 
curve). In Figs.|5]and|6]we showed the observed color-magnitude 
behavior of the photometric variations of Sa. The model magni- 
tudes are over-plotted and labeled with their respective accretion 
rates. The results of our modeling efforts can be summarized as 
follows: 



Fig. 7. Observed infrared photometry of Sa, with radiative trans- 
fer model SEDs over-plotted. The photometry is the same as in 
Figs. |5] and |6] The model SED curves are for accretion rates of 
10"^, 10"'', and 10"*^ Moyr"' (top to bottom black curve), with 
the other model parameters as given in Sect. 4.2 and Table [5] 
Also shown, with a grey curve, is the 10"^ Moyr" -model with- 
out foreground extinction. 



instantaneous variations in the temperature of the circumstellar 
material and to corresponding changes in the infrared emission 
because the dust grains are small and attain their new equi- 
librium temperature practically instantly if the radiation field 
changes, and the light travel time to the outer disk edge is only 
about 30 minutes. 

For the luminosity of the central star itself, we adopt a value 
of 10 Lq. This is lower by a factor of a few compared to the 
estimate by [Duchene et al. ( 2005 [ l, who base their value on a 
mass estimate of 2.5-3 Mq for Sa. However, dynamical model- 
ing of the Sa-Sb orbit including more recent data yields a mass 
estimate of -2.2 Mq for Sa ( |Kohler et aL 2008; Kohler 2008), 
which is essentially identical to the mass of TTauN (2.1 Mq, 
White &"Ghez 2001). The bolometric lum inosity of T Tau N i s 



relatively well determined to be 7.3^[ j Lq ([White & Ghez 



2001 



and assuming the stars are co-eval therefore yields a lummosity 
of <10Lq for Sa. A lower stellar luminosity is favorable for a 
scenario in which the IR variability is explained by variable ac- 
cretion: the required extra luminosity in "high states" is always 
relative to the stellar photosphere and thus directly scales with 
the stellar luminosity. The stellar parameters that we chose cor- 
respond to a 2 Myr-old star of 2.2 M© according to the pre-main 
sequence evolutionary tracks by Dott er et al.[ ([2008). 

We performed model calculations for mass accretion rates 
of 10"^, lO"*", and 10"^ Moyr"'. For the lowest accretion rate, 
the accretion luminosity amounts to <2% of the stellar lumi- 
nosity and the resulting fluxes are essentially identical to those 
obtained with a zero accretion rate. We calculated models for a 
range of nearly edge-on inclinations. Because our objective is to 
demonstrate that a variable accretion model is in principle capa- 
ble of reproducing the observed color-magnitude behavior, we 
searched for a model that gives a fair match to the observations 
but did not comprehensively explore parameter space. A model 
with parameters as listed earlier in this section and in Table |5] 



• the near-infrared color and brightness variations are quali- 
tatively reproduced: we obtain a "redder when fainter" be- 
havior in a K-L' vs. K color-magnitude plot (see Fig.|5]). The 
slope of the relation in the model is somewhat shallower (i.e. 
closer to grey) than observed. 

• the observed K-[12.8] vs. K color-magnitude behavior is 
qualitatively reproduced. The slope of the modeled relation 
fits the observations well. The considerable observed scatter 
can at least in part be attributed to the time lag between the 
near- and mid-infrared observations (see Fig.|6]l. 

• the disk model always produces a 10 fim silicate in emission, 
even at 90° (edge-on) inclination. Additional foreground ex- 
tinction is needed to explain the observed silicate absorption 
feature. 

Why does our model reproduce the observed "bluer when 
brighter" behavior, whereas the previous model considered in 
the context of TTauS (Ca lvet et aL][T997) [Beck et al.|[2004l i 
did not? In the Calvet model the central star is not highly ex- 
tincted. At low accretion rates, the very blue stellar photospheric 
emission dominates the near-infrared spectrum. As the accretion 
rate is increased, the comparatively red disk spectrum becomes 
more prominent and dominates the near-infrared emission above 
10"' Mgyr"'. The changing star-disk contrast causes a "redder 
when brighter" color-magnitude behavior with varying mass ac- 
cretion rate. In contrast, in our highly reddened model of the Sa 
star and disk system, the photosphere never substantially con- 
tributes directly to the near-infrared flux, because it is obscured 
by the Sa disk. At low accretion rates, the emission is domi- 
nated by the far side of the inner edge of the passive disk part, 
irradiated by the central star. With increasing accretion rate, the 
inner edge of the passive disk becomes hotter and, in addition, 
the comparatively blue active Shakura-Sunyaev disk contributes 
stronger to the near-infrared spectrum. The combined eff'ect is a 
"bluer when brighter" color-magnitude behavior. 

The variable accretion mechanism has been proposed previ- 
ously, based on the apparent grey behavior of the observed near- 
to mid-infrared flux variations ( [Ghez et al.|1991[ [Duchene et al.[ 
12002 1. However, based on a much larger data set we show that 
the large, long-term flux variations at 2.2 and 12.8 jum are ac- 
companied by substantial color variations (Fig.|6]), which can be 
reproduced by our variable accretion model. However, these data 
can also be explained with variable foreground extinction (see 
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the reddening vectors drawn in Figs. |5] and |6]). The total ampli- 
tude of the foreground extinction required to cover the observed 
brightness range is AAv~40mag. 

4.2.1. Variable extinction causing variable disk-illumination? 

In the previous paragraphs we showed that variable extinction 
in its simplest form, i.e. absorbing material passing in front of 
a source of approximately constant intrinsic luminosity, cannot 
cause the observed fast mid-lR brightness variations. However, 
the color-magnitude behavior of the large, multi-magnitude vari- 
ations that have been observed during the last two decades, can 
be explained with variable extinction (see Sect. 4.2 1. The vari- 
able accretion model is in principle capable of explaining both 
the detected sub-magnitude fast variations, as well as the pre- 
sumably slower multi-magnitude variations. Here we will con- 
sider a conceptually different realization of the variable extinc- 
tion scheme, in which we attempt to explain both the short- and 
long-term variations with extinction alone. 

To this purpose, we will introduce "central extinction 
events": variable absorption columns close to the central star 
that cause substantial changes in the illumination of the disk sur- 
face at larger radii with optical and near-infrared radiation. This 
could in principle lead to fast variations at mid-lR wavelengths 
because the material causing the variable disk illumination re- 
sides much closer to the central star than the material emitting 
the bulk of the mid-lR radiation, thereby circumventing the dy- 
namical timescale problem encountered in the "classical" vari- 
able extinction scenario. 

In order to explore this scheme in a more quantitative sense, 
we constructed a passive disk model for Sa, in which we intro- 
duce variable absorption columns by increasing the scale height 
at the inner disk edge. The disk parameters^ are the same as 
those of the model presented in Sect. |4.2| but as the central il- 
luminating source we replaced the lOL© star plus variable ac- 
cretion component with a star-only configuration with a constant 
luminosity of 40 Lq ( Duchene et al.|2005] l. 

We considered both a scenario in which the scale height at 
the inner disk edge is time variable, and one in which a fixed 
"structure" rotates at the inner disk edge with the Keplerian 
speed, casting a shadow over the outer disk that rotates along 
with the same angular velocity. For the former experiment, we 
artificially increased the scale height of the disk between and 
2Ri„, at all azimuthal angles. If we increase the scale height by 
55%, with respect to the disk parameters given in the previous 
paragraph, the flux at 12.8 fim decreases by 0.26 mag because 
the inner disk somewhat shadows the regions at larger radii. The 
aiO.26 mag brightening at 12.8 fim that we observed in February 
2008 may conversely be explained by the opposite change in in- 
ner rim structure: a reduction of the pressure scale height from 
«;0.034 AU to a!0.022 AU within <4days. The height above the 
midplane where the disk becomes optically thin to the stellar ir- 
radiatiorj^ then decreases from a!0.16AU to a;0.11AU, i.e. by 
more than twice the pressure scale height corresponding to hy- 
drostatic equilibrium. The orbital period at 0.5 AU from a 2.2 Mq 
star is a!87 days, and thus vertical movements on the order of one 
scale height require roughly 20 days, which is much longer than 
our observational constraint. Thus we conclude that changes in 



the inner disk structure cannot cause the observed fast mid-IR 
brightness variations. 

A somewhat different approach is to take a fixed enhance- 
ment of the disk scale-height at the inner rim over only part of 
the azimuthal range, causing a "shadow" to be cast over the disk 
behind it. This shadow will move over the outer disk with an 
angular velocity corresponding to the Kepler speed at the inner 
edge, and thus cause a "cool region" at larger radii that moves 
at locally vastly super-keplerian velocities over the disk surface. 
For near edge-on inclinations, where the mid-IR flux we see is 
dominated by the far side of the disk, this can cause substan- 
tial dimming. We used the 3D version of radmc (Dullemond et 
al. in prep.) to calculate light curves of a system in which the 
scale height at the inner rim is increased with respect to hydro- 
static equilibrium over a limited azimuthal rangePl With this we 
can achieve brightness variations of <0. 15 mag on timescales of 
a!20 days. Because the enhancement of the scale height is unreal- 
istically large and yet the brightness variations are substantially 
smaller and slower than observed, we do not consider this sce- 
nario to be a viable explanation for our observations. 

We conclude that variable accretion is the prime mechanism 
responsible for the observed short-term mid-IR variability, and 
exclude a significant role for variable extinction in this matter, as 
explained in Sect. |4. 1| The short-term near-IR photometric vari- 
ations are naturally also explained by variable accretion, though 
the existing data do not strictly exclude variable extinction as 
their cause. Based on the hitherto available data, both variable 
accretion and variable extinction are viable mechanisms to ex- 
plain the long-term photometric variations and the associated 
color-magnitude behavior. Data presented in past studies, in par- 
ticular a time-variable optical depth in the 3.05 yum ice absorp- 
tion feature ( Beck et al. 2001 1, do suggest that at least part of the 
variability can be attributed to variable extinction. 

4.2.2. The 10 /im silicate feature in Sa 

The fact that T Tau S shows a deep silicate absorption feature in 
the 10 yum spectral region ( Ghez et"aL]|1991[l, w hich was later 
shown to arise in Sa only ( Skemer et al.| 2008 Ratzka et al. 
2009), has been regarded as evidence that the 10 fim emission is 
dominated by a disk that is seen nearly edge-on from our vantage 
point (e.g. jDuchene et al.|2005 1. Our radiative transfer model of 
the Sa disk, however, always shows a silicate feature in emission, 
even for an edge-on inclination (see the unreddened model curve 
in Fig.|7] grey line). This is due to the "artificial" truncation of 
the Sa disk by the tidal influence of Sb, with an outer radius of 
5 AU in our model. In "normal " T Tauri disks, which are much 
larger, cold material in the outer disk regions obscures the warm 
inner disk when viewed at near edge-on inclinations, causing a 
silicate feature in absorption. In the special case of Sa, the entire 
disk is warm enough to emit in the 10 fim spectral region. The 
radiation that we receive is dominated by the optically thin up- 
per disk layers, yielding a silicate feature in emission, regardless 



of the viewing angle. Ratzka et al. (2009 1 reach the same con 



elusion based on their radiative transfer model for Sa using the 
MC3D code. Note that the shape and peak to continuum ratio 
of the silicate band depend on the grain size used. The small, 
sub-micron sized grains used in our model cause a relatively 
strong emission band peaking just short-ward of 10 /vm. The use 



7 = 2/7, 



^ i?in = 0.5AU, i?„„, = 5AU, 
tfo = 0.021 AU at 0.5 AU. 

^ Here this is taken to be the height where the optical depth in the 
radial direction atA = 0.55 fim falls below unity. 



' Over n radians in azimuth the scale height is increased with a sinu- 
soidal shape, so that at the central azimuthal angle the scale height is 
twice the hydrostatic value, and at angles of -n/2 and -t-n/2 with respect 
to the maximum it reaches the hydrostatic value. 
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of larger grains results in a weaker, broader silicate feature, and 
one can make the silicate band essentially disappear from the 
spectrum by using very large (>10 //m) grains. However, one 
never gets a silicate feature in absorption. Thus, an additional, 
external absorption component needs to be invoked. This is con- 
sistent with the fact that Sb, whose disk is thought to be seen 
more or less face-on (e.g. Duchene et al.|[2005] l, suffers an ex- 
tinction of Av~15mag. The shape of the observed silicate ab- 
sorption feature implies that the material causing the absorption 
consists predominantly of small (<1 /im) grains. 

4.3. A tentative scenario for TTauSa 

The light curves shown in Fig. |4] can qualitatively be described 
as follow^ From the beginning of the IR observations, in 1982, 
until 1989 T Tau S appears to have been in a "low state". Because 
there are only three measurements in this period, it is by no 
means certain that it has been faint during this whole period, 
but the available measurements all show TTau S to be of essen- 
tially the same (low) brightness at 2.2 jim. Sometime in 1989, 
the brightness of the southern component started to increase 
strongly. Based on data from later epochs this increase, or at least 
the vast majority of it, is attributed to Sa only. During the period 
from a;1991 to a;1995 the temporal sampling was very poor, but 
the few available measurements suggest S a to be relatively faint. 
From ^1996 to late 2001 Sa has been in a "high and variable" 
state. This period has been very well sampled by |Beck et aT] 
( [200 4') and the K-band brightness varied between Ksi7.5mag 
and Ka;6 mag. From late 2001 to late 2003 the brightness of Sa 
decreased to Ks;10.5mag. The temporal coverage of the obser- 
vations is poor and thus we cannot assess how "smooth" this 
decrease was. Sometime between late 2004 and late 2006 the 
brightness of Sa rose again to Ka!9 mag, and is currently varying 
around this value. 

In view of the variable accretion scenario, the dramatic 
brightening events in the early and late 1990s and the inter- 
spersed episodes of irregular flux variations, of which we may 
be seeing the last twitch at the current time, could be interpreted 
as an "outburst" of accretion. In this respect, it is particularly in- 
teresting to compare the timing of the outburst with that of the 
periastron passage of the companion Sb. In the current best-fit 
solution of the Sa-Sb orbit, the closest approach occurred in June 
1995, with an uncertainty of somewhat more than a year (Icr), 
as we have indicated in Fig. |4] with a grey vertical band. Thus, 
we have a period of enhanced activity around the time of closest 
approach of the companion, with the middle of the active period 
timed somewhat after periastron passage, and the strongest ac- 
tivity occurring in the second half of the active period. We will 
now argue that this behavior provides a good qualitative match 
to a scenario in which an accretion outburst in the center of the 
Sa disk is induced by the gravitational influence of Sb, tidally 
disturbing the Sa disk during periastron passage. 

When a young star and disk system has a low-mass compan- 
ion, the latter will tidally truncate the disk around the primary 
to ~l/3 times the semi-major axis of the orbit (e.g. |Artymowicz| 
1993 ). If the companion is on an eccentric orbit, its gravitational 



influence on the primary disk will be strongly time-dependent, 
being relatively weak during most of the orbit but culminat- 
ing near periastron passage. Kley & Nelson (2008 1 modeled the 
close binary system y Cephei in its youth, which is thought to 



have been very similar to the Sa-Sb system today. It has pri- 
mary and secondary masses of «1.59 M0 and a!0.38 Mq, respec- 
tively, a semi-major axis of a!l8 AU, and an eccentricity of ^0.36 
(Ha tzes et al.||2003] l. In addition, the primary has a planetary 
companion with a mass of M sin / a; 1 .7 Mjup ( Hatzes et al. 2003| l, 
which is, however, not of central importance to the current dis- 
cussion. Kley & Nelson ("2008) performed hydro-dynamical sim- 
ulations of the disk around the young y Cep primary and consid- 
ered the tidal effects of the low-mass secondary. The simulations 
reveal a nearly undisturbed disk at most orbital phases, when the 
companion is at relatively large distances. As the companion ap- 
proaches the primary, however, two strong spiral shock waves 
are excited which propagate all the way to the inner disk (see 
Fig.|8]l. This takes some time, and therefore the gravitational dis- 
turbance of the inner disk is not symmetrical in time with respect 
to the periastron passage, but lags somewhat. The strongest dis- 
turbances at the inner disk edge occur substantially after closest 
approach. 

We speculate that such shock waves induce enhanced accre- 
tion in the innermost disk region of Sa when Sb is near peri- 
astron. For example, the mechanical energy dissipated in the 
shocks may cause a normally essentially neutral disk medium 
to be sufficiently ionized for MRI to become effective, naturally 
leading to strongly enhanced accretion. We have re-run the Kley 
& Nelson model including viscous heating and radiative cooling 
for the parameters of the Sa-Sb system, and find the same behav- 
ior as seen in the yCep calculation^ The simulations provide a 
good qualitative match to the observations: a period of enhanced 
accretion activity near every periastron passage, whose peak has 
a substantial time lag with respect to the closest approach of the 
secondary. Quantitatively, the model shows much more modest 
variations in the accretion rate than inferred by the observations, 
and the time between periastron passage and maximum accre- 
tion activity is longer in the simulations by a factor of ~3. The 
Kley & Nelson model focuses on the dynamical evolution and 
mass loss of the disk at its outer edge, but does not attempt to 
treat the accretion physics at the inner disk edge in detail. A 
model with a more detailed treatment of the disk viscosity and 
accretion physics in the innermost disk regions is needed, which 
allows the spiral shock waves to propagate inward more quickly 
and can account for the required high accretion rates. This is be- 
yond the scope of the current work, and whether our qualitative 
model can be refined to provide a quantitative match to the ob- 
servations once all the relevant physics are included remains to 
be confirmed. 

Here we pose the companion-induced accretion outburst as 
an idea to explain the observed behavior of T Tau S during the 
past decades, and provide a clear prediction that can be easily 
verified observationally: if this scenario is correct, Sa should re- 
turn to its "low state" some time in the near future, and stay 
there until the next periastron passage of Sb. Exactly how long 
it will take for Sa to return to "quiescence" in detail depends on 
the actual accretion physics and inner disk parameters, neither 
of which we know. Thus, while it is technically trivial to test the 
proposed scenario with currently available observing facilities, 
it may in practice take a substantial number of years to convinc- 
ingly falsify our scheme. Supportive evidence could be built up 
over the years if Sa is observed to remain faint, the final prove 
may come only many years from now if the next outburst co- 
incides with the next periastron passage. The current astromet- 



We will describe the K-band light curve because it has the best tem- 
poral sampling. The L' and 12.8 fim light curves follow the K-band (see 
also Figs.|5]and[6]l. 



' Note that the model discussed here and shown in Fig.jsjis a hydro- 
dynamical calculation and is independent of the radiative transfer model 
discussed in Sect. 4.2 



14 



R. van Boekel et al.: Variable accretion as a mechanism for brightness variations in T Tau S 




Fig. 8. Two-dimensional temperature distribution of the Sa disk 
10.8 yrs after the periastron passage of Sb, for a hydrodynamic 
simulation using the physical parameters of the T Tau S system. 
X and Y are spatial coordinates, and the disk is viewed face- 
on in this representation. To better bring out the entire range of 
temperatures we show T[K]^^'^, the colors cover the temperature 
range of 10 to 2000 K. 



ric data still allow a wide range of orbital periods, from a!25 to 
~100 years, but as new observations accumulate over the com- 
ing years, we should be able to fairly accurately predict the onset 
of the next outburst long before it would occur. 

4.4. Issues 

Here we briefly discuss some issues that are relevant to the cur- 
rent discussion. 



4.4.1 . The mass and lifetime of tine Sa disk 

The amount of material that we can reasonable assume to be 
present in the Sa disk is at most a few times 10"^ Mq. If the 
average accretion rate during an outburst is 5 x lO"*" Moyr"' and 
an outburst lasts for 10 years, we can have only on the order 
of lO'^ tiiese events before the disk is all gone, unless the Sa 
disk can be supplied with fresh material. If the orbital period 
of Sa-Sb is 50 yrs, this means the current situation can exist for 
only <5000yrs. Possibly the Sa-Sb pair became a close binary 
only within the last several thousand years due to gravitational 
interaction with a third object. It may also be possible that the 
Sa disk is still supplied with fresh material coming in at high 
latitudes with respect to the orbital plane as the Sa-Sb pair moves 
through ambient cloud material. Because the disks of Sa and Sb 



appear strongly misaligned (Duchene et al. 2005 Ratzka et al 



4.4.2. A note on tine Sb disk 

Like the primary star Sa, the low-mass companion Sb also shows 
a substantial infrared excess, indicating that some circumstellar 
material is still present, presumably also in a small circumstel- 
lar disk. As our calculations show, Sb will strongly disturb the 
Sa disk during periastron passage. Conversely, the gravitational 
pull of Sa will affect the disk around Sb even more strongly, 
and accretion onto Sb should also be enhanced during periastron 
passage. Therefore one may wonder why the observed bright- 
ness variations are much smaller in Sb (maximum amplitude 
aiO.4 mag RMS -0.2 mag) than in Sa (amplitude >3 mag). 

|Duchene et al.| ( |2002 2005 1 ) measured the Bry line flux 
of Sb to be 7.2xir^in=^on 2000 November 19, and 
9.4x10"'^ Wm"^ on 2003 December 12, with corresponding K- 
band continuum magnitudes of 8.77 and 8.50, respectively, the 
spectral type of Sb corresponds to reff~3800K (Duchene et al. 



2009[ l, the Sa disk and the Sa-Sb orbit are misaligned, and cloud 
material is present close to the Sa-Sb system (Se ct.|3 .2.3||Herbst| 
let al.|2007HBeck et al.|2008l|Gustafsson et al.|2008) , the notion 
that cloud material is still being freshly supplied to the Sa disk 
is realistic. 



2005 ), and assuming an age of 2Myr this yields a mass of 
approximately 0.6 M© and a luminosity of a;0.45Lo by com- 
parison to the pre-main sequence evolutionary tracks of Dotter 
et al. ( 2008 1, which seem reasonable. The corresponding stel- 
lar radius of Sb is /?,=1.55R0. Assuming a foreground extinc- 
tion of Av=15mag, a standard extinction law (Ak=0.1 12 Ay), 
the relation between Bry luminosity and accretion luminosity 
given in equation 3 of |Calvet et al.| ( |2004j ), and using the ap- 
proximate relation Lacc=GM*Macc/^*, one derives mass accre- 
tion rates of 3.5x10"^ Moyr"' and 4.4x10"^ Moyr"' for the re- 
spective epochs. The corresponding accretion luminosities are 
approximately equal to the stellar luminosity. Therefore, variable 
accretion is a viable mechanism for the observed modest near- 
infrared variations in Sb as well. How variations in the accretion 
luminosity translate to variations in the near-infrared continuum 
excess emission is somewhat dependent on the disk structure: the 
extra luminosity is released in the UV which cannot be observed 
due to the high extinction, and we only see the reprocessed radi- 
ation from the hot inner disk regions. 

We note, however, that the existing data also perfectly agree 
with a time variable extinction towards Sb. Between the afore- 
mentioned epochs of Bry spectroscopy by Duchene et al. the 
extinction would have decreased by AAva!2.4mag, assuming a 
standard extinction law. 



4.4.3. Mid-infrared briglitness during tine peak of tine outburst 

We are not aware of mid-lR measurements performe between 
roughly 1998 January and 2001 November, but we argue Sa 
must have been very bright during that period, peaking at »40 Jy 
at 12.8 yum. Should any mid-lR observations have been per- 
formed during this period, they are highly relevant to the cur- 
rent discussion. The variable extinction scenario predicts a lower 
peak-brightness in mid-lR than the variable accretion scenario. 
However, only if mid-lR observations have been performed 
(near-) simultaneously with the existing K-band data they may 
provide supportive evidenced for either scenario. In case of a 
substantial time lag the potential short-term variations will in- 
hibit drawing firm conclusions. 



5. Summary 

The southern infrared companion of T Tau is a close binary sys- 
tem of which the more massive member, TTauSa, shows strong 
long-term photometric variability, with an amplitude of >3 mag 
at 2.2 fim and a!2mag in the 10 fim region on timescales of 
years. Short-term variability of ailmag within one week has 
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been detected at 2.2 /im. The physical mechanism driving the 
observed flux variations has been under debate, the originally 
proposed intrinsic luminosity variations due to a variable accre- 
tion rate have later been challenged in favor of apparent bright- 
ness fluctuations due to variable foreground extinction. 

We revisited the nature of the photometric variability of 
T Tau Sa using new observations in the mid-infrared. In February 
2008 we detected a very fast increase of the mid-IR brightness 
of TTauS which can be attributed to Sa: +26+2% at 12.8/imin 
3.94 days. Using simple geometric arguments and basic physic 
laws we can exclude that this short-term variability is due to 
time-variable extinction. We show that the fast brightness fluc- 
tuations must instead reflect changes in intrinsic luminosity, and 
attribute these to variable accretion. The key aspect of our argu- 
ment is the short timescale of the variations in combination with 
the long observing wavelength. 

We show that variable accretion can plausibly explain also 
the much larger long-term variability and its associated near- 
infrared "bluer when brighter" color-magnitude behavior, by 
creating a 2D radiative transfer disk model of the Sa disk that 
qualitatively reproduces the observed properties. However, the 
hitherto observed long-term variability can also be explained 
within the variable foreground-extinction scenario. A combina- 
tion of both mechanisms, variable accretion and variable fore- 
ground extinction, may be required to explain the collective pho- 
tometric variability of T Tau Sa. 

In view of the tentative scenario that both the short- and long- 
term variability of Sa are caused by variable accretion, we re- 
viewed the infrared light curve in relation to the parameters of 
the Sa-Sb orbit. We propose that the dramatic brightening events 
in the early and late 1990s and the interspersed episodes of irreg- 
ular flux variations were an accretion outburst, induced by gravi- 
tational perturbation of the Sa disk during the periastron passage 
of Sb, with closest approach in »il995. The induced maximiun 
accretion rate is on the order of 10"^ Moyr"'. While we do not 
model the physics of this outburst in detail, we make a clear ob- 
servational prediction: if this scheme is correct, Sa should return 
to "quiescence" some time in the near future, and remain in a 
low state until the next periastron passage of Sb. 
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